Dynamics and thermodynamics of fragment emission from excited sources 
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In this paper we study the process of fragmentation of highly excited Lennard- Jones drops by 
introducing the concept of emitted fragments (clusters recognizable in configuration space which 
live more than a minimum lifetime). We focus on the dynamics and thermodynamics of the emitted 
sources, and show, among other things, that this kind of process can not be cast into a sequential or 
simultaneous one, and how a local equilibrium scenario comes up, allowing us to define and explore 
a local temperature, which turns out to be a strongly time-dependent quantity. 
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I. INTRODUCTION 

The problem of fragmentation attracts the interest of 
physicists in many branches of physics. In particular, in 
this problem we have to deal with systems that are both 
finite and non-extensive. Moreover, we can face from 
extremely small systems like the nucleus to huge ones 
like galaxies. In both cases the size of the system under 
analysis is of the order of the range of the interaction 
potential and calls for a careful analysis. In what follows 
we will focus on the analysis of small Lennard- Jones [LJ) 
drops i.e. a system in which the dynamics is generated 
by a Hamiltonian in which there is a very short-range 
strong repulsion and a short-range attraction. As such we 
consider that it will give relevant insight into the problem 
of fragmentation. 

This problem has been analyzed from different points 
of view and we classify them in two main groups i.e. Sta- 
tistical models and dynamical models. In the first case 
Q 01^ the main assumption is that the fragmentation 
process is driven by phase-space occupancies. The always 
present assumption of freeze-out volume states that the 
system fragments in an equilibrium scenario at a given 
fixed volume and that thermal and chemical equilibrium 
is reached, in correspondence with the flavour chosen by 
the researcher (i.e. microcanonical, canonical or grand- 
canonical). Recently this kind of analysis has been ex- 
tended to the isobaric isothermal ensemble Even in 
this kind of formulation there have been some efforts to 
include non-equilibrium effects, like for example the pres- 
ence of expansive collective motion 

The other kind of models we dubed, dynamical ones, 
include among others: Classical ones JJ, quantual 
ones 0, which are fully microscopic. In this cate- 
gory we also include those based on numerical realiza- 
tion of kinetic equations like BUU, VLD, etc. Focusing 
on the fully microscopic models the main advantage we 
find is that we have correlations of all orders at all times. 
Moreover non-equilibrium features of the process can be 
readily explored. 

Continuing with a series of previous works we will 



study the dynamics of highly excited 147 particles LJ 
drops 7]. Previously we have mainly considered the 
problem of fragment formation. We have shown that this 
phenomena takes place in phase space and as such it is 
not an experimental observable. In this communication 
we will focus our attention on the properties of the emit- 
ting sources and the corresponding emitted fragments. 
At this point it is worthwhile to define what we under- 
stand as emission. We say that a fragment has been 
emitted when it is recognizable in configuration space. 
We have already made this classification S'lid we 

have shown that there are well defined time scales for 
the processes of fragment formation and fragment emis- 
sion. 

In Section // we review the different fragment recog- 
nition algorithms currently in use. In Section /// we 
present different Temperature definitions that will be rel- 
evant in the analysis of the evolution of the excited drops. 
Section IV deals with the model we study. Section V in- 
cludes all the results we have obtained in our numerical 
simulations (fragment mass distributions, characteristic 
times, temperature of the emitting sources, role of the 
radial flux, etc) . Finally conclusions will be drawn. 



II. FRAGMENT RECOGNITION ALGORITHMS 

The problem of analyzing molecular dynamics calcu- 
lations is an old one and is not completely settled. To 
our knowledge there are three main fragment recognition 
algorithms in use: MST, MSTE, and ECRA. 

The simplest and more intuitive cluster definition is 
based on correlations in configuration space: a particle i 
belongs to a cluster C if there is another particle j that 
belongs to C and |ri — rj| < rd, where rd is a parameter 
called the clusterization radius. If the interaction poten- 
tial has a cut off radius Vcuti then rd must be equal or 
smaller than Vcut, in this work we chose rd — Vcut = 3(7. 
The algorithm that recognizes these clusters is known as 
the "Minimum Spanning Tree" (MST). The main draw- 
backs of this method is that only correlations in q-space 
are used, neglecting completely the effect of momentum. 
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An extension of the MST is the "Minimum Spanning 
Tree in Energy space" (MSTE) algorithm iljj. In this 
case, a given a set of particles i,j,...,k, belongs to the 
same cluster Ci if: 

VieC, , 3jeC,/e,j <0 (1) 

where = V{rij) + (p^ — pj)^/2/x, and ^ is the re- 
duced mass of the pair MSTE searches for config- 
urational correlations between particles considering the 
relative momenta of particle pairs. In spite of not be- 
ing supported by a physically-sound definition of a clus- 
ter, the MSTE algorithm typically recognizes fragments 
earlier than MST. Furthermore, due to its sensitivity in 
recognizing promptly emitted particles, it can be useful 
to study the pre-equilibrium energy distribution of the 
participant particles. 

A more robust algorithm is based on the "Most Bound 
Partition" (MBP) of the system The MBP is the 

set of clusters {Ci} for which the sum of the fragment 
internal energies attains its minimum value: 

{C^} = argmm[E^c,} ^J2^EZ't] 
{c\} 

= E[E^^r+ E ^M-] (2) 

j<k 

where the first sum in is over the clusters of the par- 
tition, is the kinetic energy of particle j measured 
in the center of mass frame of the cluster which contains 
particle j, and Vij stands for the inter-particle poten- 
tial. It can be shown that clusters belonging to the MBP 
are related to the most-bound density fluctuation in r-p 
space 0. 

The algorithm that finds the MBP is known as the 
"Early Cluster Recognition Algorithm" (ECRA). Since 
ECRA searches for the most-bound density fluctuations 
in q-p space, valuable space and velocity correlations can 
be extracted at all times, specially at the very early stages 
of the evolution. This has been used extensively in many 
fragmentation studies 

II Sill 13 and has helped to 
discover that excited drops break very early in the evo- 
lution. 

When we use the three above-mentioned algorithms 
and we apply a criterium based on the average micro- 
scopic stability of the clusters (see for example y) to 
determine the corresponding times of fragment forma- 
tion three time scales emerge, which satisfy the following 
relation: tecra < tmste < tmst- But the meaning 
of each of these times is quite different, tecra refers 
to that time at which, on average, clusters attain mi- 
croscopic stability regardless of the structure of the sys- 
tem in q ~ space, tmste has a rather obscure meaning 
because MSTE algorithm is not well defined from the 
physical point of view for dynamical problems. Finally 
Tmst refers to that time at which, on average, micro- 
scopic stability is attained for free fragments. Please take 
into account that depending on the definition of rd this 



last condition can be rephrased as "weakly interacting 
fragments" . 

III. TEMPERATURES 

Because we are aiming to perform some kind of ther- 
modynamic analysis of the system, a prescription for the 
calculation of the temperature in small (and probably out 
of equilibrium) systems, is to be given. We should keep 
in mind that we are going to analyze physical processes 
in the framework of the molecular dynamics ensemble, 
which is the microcanonical ensemble plus the conserva- 
tion of total momentum. Let us assume that at a given 
stage of the evolution of the system we can identify the 
biggest fragment according to MST algorithm. 

In order to study the time evolution of the thermody- 
namics of the biggest emiting source we define the local 
temperature Tioc{t) as the velocity fluctuations around 
the mean radial expansion of the biggest source. So we 
first introduce the mean radial velocity of the biggest 
source, which is defined as 

Where Nbf indicates the number of particles of the 
biggest source at time t, and <>e denotes an average 
over all events. Both Vi(t) and Vi{t) are measured from 
the center of mass of the biggest source. Taking into ac- 
count that collective motion should not be considered to 
calculate the temperature, we define the local tempera- 
ture as Tioc'- 

2 m 

Tlocit) = E - ''radit) ■ Ut)f (4) 

2—1 

This definition is closely related to the one used in 
for the analysis of expanding systems. In that case 
we obtained the local temperature of the system by the 
following procedure: We divided our drops in concentric 
spherical regions, centered in the cm of the system, of 
width 5r = 2a. The mean radial velocity of region i in 
this case is: 

where the first sum runs over the different events for a 
given energy, the second over the particles j that belong, 
at time to region i; Vj and are the velocity and 
position of particle j. iVi(f) is total number of particles 
belonging to region i in all the events. Then the local 
temperature, which will be called T^^^^u, is defined as: 

(.) _2 1y^l / ^iVl-.Y (6) 
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where Ni is the total number of particles in cell i in all 
events. 

The validity of Eq^ and EqEl rely on the conjecture 
that the fragmenting system achieves local equilibrium 
(local equilibrium hypothesis). Whereas in the last case 
all particles in the inner shells should be equilibrated, the 
advantage of dealing with the local temperature (Eq^J 
relies on the fact that only the biggest source is assumed 
to have reached some degree of equilibration. In this pa- 
per we will focus mainly on this definition. If the biggest 
source is assumed to have reached equilibrium, the veloc- 
ity distribution should follow: 

/(-)=p(f^)^/^e-''¥(— )^ (7) 

ZTT 

A first approach on the study of the local equilibrium 
hypothesis (LEH) consists in analyzing the isotropy of 
the velocity fluctuations around the expansion. For this 
purpose we introduce the radial and transversal local 
temperatures: 

Nbf 

2—1 

A more profound analysis of the accuracy of the LEH 
is the study of the velocity distribution function and its 
comparison with (Eq[7I). In order to perform this study 
we are to check the significance S of the fit of the distri- 
bution numerically obtained with the EqEI To address 
the question if both distributions are different we will 
perform perform a Pearson test ff^. If the signifi- 
cance is high enough not to reject the hypothesis that 
both distributions differ we obtain as a by product the 
possibility of calculating the temperature as: 

Tmax = ma^ (10) 
Where a is the width of the velocity distribution. 

IV. THE MODEL 

The system under study is composed, as in previous 
works, by = 147 particles interacting via a truncated 
and shifted Lennard-Jones (LJ) potential, with a cut off 
radius Vcut = Scr. Energies are measured in units of the 
potential well (e), and the distance at which the LJ po- 
tential changes sign (cr), respectively. The unit of time 
used is to = y^ahn/ASe. 

The typical frequency of such potential is I'o ^ O-^'^ • 
This defines a minimum time scale for the stability of an 
interacting system. 



The equations of motion were integrated using the ve- 
locity Verlet algorithm, which preserves volume in phase 
space. We used an integration time step of O.Olto, and 
performed a microcanonical sampling every Uq up to a 
final time of t = 250io. 

This time scale was chosen because is has proven to 
be long enough as to the system to attain microscopic 
stability. 

As stated in the previous section, we are mainly in- 
terested in the properties of fragments in configuration 
space. Then we are to analyze the AID evolutions using 
the MST algorithm. We define a fragmentation process 
when a source emits a stable fragment of at least 4 parti- 
cles, i.e. the MST determination of fragments is comple- 
mented with a temporal stability condition. Otherwise 
we will be facing an evaporation process. The evolutions 
are analyzed in an event-by-event basis and the times at 
which fragmentation takes place are determined. 

V. NUMERICAL EXPERIMENTS 

Initial configurations are built as dense drops in q — 
space (region C of the phase diagram ^^). The total 
linear and angular momentum of the drops are removed 
and then, velocities are rescaled with a Maxwellian distri- 
bution so that the system has any desired value of total 
energy. 

We covered a broad range of energies, namely — 2.0e, 
-l.Oe, -0.5e, -0.2e, O.Oe, 0.2e, 0.5e, l.Oe and 2.0e (all 
reported energies are per particle). For each energy 1000 
events were calculated. 



A. Fragment mass distributions 

The range of energies we used covers all regions of in- 
terest in multifragmentation. In Fig. ^ it can be clearly 
seen that the asymptotic mass distributions go from " U- 
shaped" , at low energies to exponentially decaying ones, 
at high energies. Between these shapes a power-law 
(panel b) in Fig. ^ can be found. Technical details of 
the fitting procedure can be found elsewhere ^3 but we 
would like to point out that the fitting procedure excludes 
the biggest fragment at each event, and also explicity ex- 
cludes monomers, dimers and trimers. 

We have found interesting to study the multiplicity 
distribution of the number of times the system fragments 
for all events as a function of the energy of the system. 
In this case our approach is to follow the dynamics of the 
biggest source. At a given time ti — nto we identify the 
biggest MST fragment and check for the particles that 
belonged to this fragment at ti, and formed a smaller 
cluster at t — {n + 5)to- We have chosen 5to as our 
threshold to accept the cluster as stable because that is 
of the order of the inverse of the natural frequency. In 
this way, we show in Fig[21 this distribution for energies 
per particle E = — 2.0e, E — O.Oe, and E = 2.0e. It can be 
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FIG. 1: Asymptotic mass spectra. Energies 
E = -2.0e,E = O.Oe, and E = 2.0e are 
displayed in panels a), b), and c), respectively. 



FIG. 2: Multiplicity distribution for E = 
-2.0e,£; = O.Oe, and E = 2.0e are 
shown in panels a), 6), and c), respectively. 



easily seen that, in this range of energy, as the energy of 
the system is increased, the maximum of the multiplicity 
distribution shifts towards higher values while the width 
of the distribution increases. 



B. Characteristic times 

Some discussion regarding the nature of the fragmenta- 
tion process has focused on the distinction between simul- 
taneous and sequential kind of phenomena. In our cal- 
culation the analysis of these kind of things are straight- 
forward. Anyway, there are different characteristic times 
which can be explored. First we will look at the first 
and last time of emission. This means that we record for 
each of our evolutions the time at which the first frag- 
mentation takes place according to the aboved- mentioned 
definition of an emission process. The result of that cal- 
culation are displayed in Fig|31 

We can see that at high energies {E — 2.0e) the distri- 
butions of times are narrower that at lower energies and 
does not seem to fit the view of a simultaneous process in 
configuration space (This issue will be further analyzed 
in Section lV G|l . Another view of the same data is shown 
in Fig^in which what we show is the corresponding lapse 
of time between the first and last emission for each event. 
Once again, the distribution gets broader as the energy 
gets lower (the comparison should be restricted to events 
corresponding to E = O.Oe and E = 2.0e because of the 
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FIG. 3: (Color online) Distribution of first (full 
line) and last (dotted line) time of emission. En- 
ergies: E ^ -2.0e,E = O.Oe, and E ^ 2.Ge 
are displayed in panels a), 6), and c), respectively. 



5 



very few emission processes that are found at E = — 2.0e. 
In order to further illustrate the kind of process we are 
facing we show in FigElthe distribution of the mass of the 
emitting sources for the times of first and last emissions. 
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FIG. 5: (Color online) Distribution of mass of the emit- 
ting sources for the first (full line) and last (dotted line) 
emissions. Energies: E = —2.0e,E = O.Oe, and E = 
2.0e are displayed in panels a), b), and c), respectively. 



FIG. 4: Distribution of lapse of time be- 
tween the first and last emissions. Energies: 
E = -2.0e,E = O.Oe, and E = 2.0e are dis- 
played in panels a), 6), and c), respectively. 

The most relevant issue of this figure is the presence, 
in panel b), of a broad range of masses at which the last 
emission takes place, since the value of energy {E = O.Oe) 
correspond to a power-law mass distribution. 

The reason why we stated above that these results do 
not seem to fit the simultaneous emission picture without 
assuring it, is because we can also explore the possibility 
of occurrence of massive emission events. Even though 
emission events can be found in a rather large scale of 
time one may wonder if massive emission events have a 
definite characteristic time which would sustain the si- 
multaneous picture. 

We have used two definitions of "massive emission 
events (MEE)". The first one considers a MEE if the 
mass of the emitted fragment is at least 30% of the total 
mass of the biggest source at the time of emission. The 
second one states that a MEE has taken place if the mass 
of the biggest emitted fragment is at least composed of 
40 particles. The results of such calculations are shown 
in Figini 

It is immediate that massive emission, according to 
definition 1, can take place at any time during the evo- 
lution, while when using the second definition (which is 
more restrictive) there are barely no massive events for 
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FIG. 6: (Color online) Distribution of emission 
times of massive fragments using criterium 1 (top) 
and criterium 2 (bottom); see text for details. 

times greater than t ^ lOOto. Also notice that the lower 
value of energy we used {E = — 2.0e) is not present be- 
cause no massive fragmentations were found. 

In conclusion, when gathering all these results together 
we reach the conclusion that for the energies displayed in 
these figures the process can be viewed as a mixture of 
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sequential and simultaneous break up. 



C. Temperatures of the emitting sources 

It is natural to think that as the sources emit frag- 
ments they will undergo a cooling process. We should 
keep in mind that we are trying to analyze things in an 
event-by-event basis, without performing averages that 
would obscure the picture. In order to gain knowledge of 
the microscopic view of the fragmentation phenomena, 
we have calculated the evolution of the temperature, for 
a given energy and for each time step in which a frag- 
mentation event takes place. The temperature is the one 
of the emitting source in the time step prior to the one 
in which the emission takes place. We then plot all this 
information in a single graph (see Fig[7I). 




FIG. 7: (Color online) Internal kinetic energy of the 
emitting source as a function of its total internal en- 
ergy for three ranks. Empty (red) squares denote 1st 
emission process. Empty (green) circles denote the 
8th emission process, and empty (blue) triangles indi- 
cate the 12th emission process. Full (black) triangles 
stand for the caloric curve of the expanding system. 



kind of measurements we are referring to at this point 
do not correspond to stable systems. What might be im- 
proper is to talk about temperature instead of "effective 
temperature" . A more quantitative calculation is under 
progress and will be comunnicated shortly. 



D. Tlie role of radial flux 

In a series of previous works we have shown that, 
if the presence of collective (nonthermal) motion is not 
taken into account one find a wrong result: The presence 
of a vapor-branch in the caloric curve. Only when the ra- 
dial flux is properly incorporated in the definition of the 
temperature one find that the nonequilibrium process of 
multifragmentation appears as an almost constant tem- 
perature region at high energies. In particular, in paper 
UJ we have calculated the time of fragment formation 
of this system, (using the ECRA phase-space method). 
The values of the mean radial velocity calculated for the 
biggest source in each event according to Eq|31are shown 
in FigEl 
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FIG. 8: Radial flux as a function of en- 
ergy at time of fragment formation t// 



In Fig[7| we show the result for the analysis of 1000 
events at an initial energy of E' = 2.0e. In order to 
make the figure more readable we are only showing the 
above-described temperature for the first, last and maxi- 
mum multiplicity fragmentation events. It is clearly seen 
that the internal kinetic energy is linear with the total 
internal energy and displays a cooling behavior. This 
overall behavior will not change if we consider not only 
fragmentation events but also evaporation events. It is 
interesting to notice that if one is to calculate temper- 
atures from the analysis of the emitted light fragments 
(monomers, dimers, trimers) one would be sampling a 
source that starting from a rather high temperature, 
cools down monotonically. So the corresponding tem- 
perature should display a maximum at early times and 
decrease later on ll8|. This is not a problem at all be- 
cause the system is out of equilibrium and as such the 



In Figl^lwe show the local temperature of the biggest 
source at time of fragment formation, as a function of 
the total energy (full line) The dotted line is included in 
order to emphazise the effect of not removing the radial 
collective motion when calculating caloric curves. As ex- 
plained in Section fun temperatures were calculated from 
velocity fluctuations around the collective motion. 

In order to see how the picture emerging from a phase 
space analysis relates to the main topic of this paper 
i.e. the analysis in configuration space, we proceeded to 
calculate the temperature of the emitting sources at the 
stage of the evolution which correspond to the maximum 
of the multiplicity distribution (for example, looking at 
FigEb) the maximum of the distribution turns out to be 
rank=4) . From this kind of analysis we get the empty cir- 
cles of FigEJl It can be easily seen that both approaches 
give the same temperature. 
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FIG. 9: (Local temperature) Local temperature 
of the biggest source at Tff (straight red line) 
and "fake temperature" (dotted blue line). The 
fake temperature is calculated from the total ki- 
netic energy of the biggest source without tak- 
ing into account the importance of the radial flux 
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FIG. 10: (Color online) In this figure we show a 
comparison of caloric curves obtained by different ther- 
mometers. Empty circles correspond to an event- 
based temperature. Filled squares correspond to the 
temperature of the biggest source at time of fra- 
ment formation. Triangles up denote the local tem- 
perature (Tsheiis) and the dotted line with down tri- 
angles corresponds to the caloric curve obtained for 
a constrained system with a density p = 0.17a^'^. 



E. Local Equilibrium Hypothesis (LEH) 

In previous sections we have made use of the hypothesis 
of local equilibrium. 

In order to check this hypothesis we have performed 
the following calculation. Because in a local equilibrium 
scenario the fluctuations of velocity should be the same 
for the radial and transversal directions, we evaluated 
Eqs. IHl and O for a single event at the three energies 
chosen as examples in this work. In each of the three 
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FIG. 11: (Color online) Radial (dotted line), transver- 
sal (thin line), and local (full line) temperatures 
of the biggest source for E = — 2.0e (panel a)), 
E = O.Oe (panel bj), and E = 2.0e (panel c)) 



panels of Fig llll each corresponding to a different total 
energy (see caption for details), we show the three defi- 
nitions of temperature above-mentioned. It is immediate 
that all of them are essentially the same, except for very 
early stages of the evolution (times much shorter than 
the time of fragment formation). This suggest that the 
the LEH at time of fragment formation seems plausible 
for all energies. 



F. Maxwellian distribution of velocities 

In this section we will show that the velocity distribu- 
tion of the biggest source at time of fragment formation is 
indeed Maxwellian. Moreover, the temperature obtained 
from the standard deviation of the velocity distribution 
is almost exactly the same of Eq01 showing the consis- 
tency of our numerical studies. In Fig ll2l we show the 
histogram of the velocity distribution at time of fragment 
formation and its Maxwellian fit for E = 2.0e. 

In adittion to the excellent agreement between the ve- 
locity distributions and their fits we performed a Pear- 
son test, trying to reject the hypothesis that the 
distribution of velocities and that corresponding to the 
Maxwellian fit differ. We found a significance (S) above 
0.25 for all cases, showing that, even with a very low 
confidence level like CL = 0.80 we could reject the hy- 
pothesis that both distributions differ (i.e. the velocity 
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FIG. 12: (Color online) Velocity distribution (his- 
togram) and Maxwellian fit (dotted line) for E = +2.0e 

distribution is indeed Maxwellian). 

TABLE I: Comparison of thermometers and 
statistical significance of the Maxwellian fit 







Tfuie) 






S 


-2.0 


0.627 ± 0.043 


0.623 


0.645 


0.647 


0.86 


0.0 


0.619 ± 0.077 


0.630 


0.608 


0.674 


0.82 


2.0 


0.587 ± 0.062 


0.575 


0.506 


0.558 


0.78 



quality of the resulting fit is remarkable indeed, specially 
when one consider that we are facing an out of equilib- 
rium non-simultaneous process while the very nature of 
the binomial process requires a constant value of p for 
the whole emission process. 
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FIG. 13: (Color online) MultipHcity distribution (his- 
tograms) and binomial fit (dotted line). For energies 
(from top to bottom and left to right): E = — 2.0e, E = 
~0.2e, E = O.Oe,E = 0.2e, E = l.Oe, and E = 2.0e 



In Table./ we show a cross comparison of different ther- 
mometers. In addition to the already presented results 
we add the temperature obtained from the Maxwellian 
fit of the velocity distribution, and we also show the ob- 
tained statistical significance. 



In order to further illustrate the accuracy of this ap- 
proach we show in Fig^] the calculated probability of 
emitting n fragments during the entire process as a func- 
tion of the total energy, whith n = — 5. P{n) is cal- 
culated asssuming a binomial distribution ( Eq. with 
the value of p obtained from the best fit. 



G. Reducibility 

Not long ago, i t has been proposed by Moretto and 
co-workers [l^, ll^ that the complex process of fragment 
emission could be described in terms of a binomial dis- 
tribution. This approach rests on the assumption that 
a single transition probability p is capable of describing 
the emission process when no regard is pay to the mass 
or composition of the emitted fragments. In this way, the 
probability of emitting n fragments in a series of m trials 
should follow the well-known binomial distribution: 
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(11) 



In which m stands for the number of "trials", while n 
stands for the number of successes. Following 19], we 
associate the parameter m to the maximum multiplicity 
for each energy and p stands for the transition probabil- 
ity. In Fig^lwe show the result of such an analysis. The 
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FIG. 14: (Color online) Probability of emit- 
ing (red), 1 (green), 2 (blue), 3 (violet), 4 
(cyan), 5 (yellow) fragments as a function of energy 

We will not further analyze the implications of this 
binomial fit because we have not been able to calculate 
microscopically the associated transitions barriers yet. 
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VI. CONCLUSIONS 

In this work we have analyzed the dynamics and ther- 
modynamics of the emission process in a classical frag- 
menting system. We have focused on the evolution of 
the biggest source, and the corresponding emitted frag- 
ments. We have shown that the emission process can not 
be cast into a sequential nor into a simultaneous scenario, 
because even though there are many events in which frag- 
ments are emitted sequentially the phenomenon of mas- 
sive emission is frequent enough as to forbid the char- 
acterization as a purely sequential one. We have also 
shown that when looking at fragments well defined in 
configuration space, the process of emission can not be 
cast into a isothermal one, i.e. temperatures are time- 
dependent. Moreover, there is a time-dependence of the 
size of the emitting source. As such, standard thermo- 
dynamical models which fix temperature and volume do 
not seem to be appropiate to reveal the true nature of 
the phenomenon under analysis. We have also shown, 
by analyzing the temperature of the source at the stage 
of most probable multiplicity, that it is possible to re- 
cover the caloric ciirve already obtained in the frame of 
phase-space analysis. If we recall Table / we see that 
temperature can be determined according to a few rea- 
sonable definitions, which provides in all cases consistent 



results. 

Two other interesting results have been obtained as 
a by-product of these calculations: First we should like 
to mention that the velocity distribution functions of 
the particles that form the biggest source are remark- 
ably Maxwellian for times larger than the time of frag- 
ment formation after the collective expanding mode is 
removed. Second we have found that the probability of 
fragment emission (summed-up over all sizes) is amaz- 
ingly binomial, suggesting a constancy of the transition 
amplitudes. We do not have a microscopic description of 
this constancy right now and we are currently working 
on it. 

We hope that these findings will encourage the devel- 
opment of new, more accurate and realistic models to 
describe nuclear multifragmentation. 
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